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Abstract 

A dual-stacked patch antenna is analyzed using a hybrid finite 
element - boundary integral (FE-BI) method. The metallic patches 
of the antenna are modeled as perfectly electric conducting (PEC) 
plates stacked on top of two different dielectric layers. The antenna 
patches may be of any shape and the lower patch is fed by a coaxial 
cable from underneath the ground plane or by an aperture coupled 
microstrip line. The ability of the hybrid FEM technique for the 
stacked patch antenna characterization will be stressed, and the EM 
coupling mechanism is also discussed with the aid of the computed 
near field patterns around the patches. 


I. Introduction 

Since being proposed, the dual stacked patch antennas received tremendous 
attention mainly because of the drastic improvement of their operational 
bandwidth in comparison to the inherent narrow band feature of single patch 
configurations. Further bandwidth improvement was also reported by em- 
ploying an aperture coupled microstrip line, where the lower patch is eiectro- 
magnetically coupled to a microstrip line through an opening at the ground 


58 


plane. This design was shown capable of achieving over 10% bandwidth. 
Certain numerical approaches for extensive analysis of non-rectangular sin- 
gle patch antennas with an aperture coupled feed have been investigated 
(e.g. (1, 2]). Unfortunately reports for non-rectangular stacked structure 
analyses have not been seen in the literature, except for experimental results 
published in [3, 4]). An analysis was only reported in [5] where the Hankel 
transform was employed on the assumption that the feed structure is absent 
from the analytical system. This implies that any change in the resonance 
due to the feed can not be modeled using this technique. A comprehensive 
simulation for stacked non-rectangular patch configurations is therefore de- 
sired for design purposes. This is the motivation of this report and we will 
discuss both probe and aperture coupled feeds. The former offers flexibility 
in the feeding network layout and makes the feed design separable from that 
of antenna elements [3], whereas the latter design, while challenging for most 
of the existing methodologies, is readily incorporated into our proposed hy- 
brid FEM formulation, especially when the patches and the coupling opening 
(aperture) are non-rectangular for wider bandwidth designs. 

In the following sections, we first briefly describe the hybrid FEM formu- 
lation as applied to the patch antenna analysis. More detailed discussions 
of the numerical approach may be found in [6] and [7]. Using our numerical 
simulator for characterizing the antenna’s dielectric and geometric parame- 
ters, we concentrate on the microstrip feed line modeling with the aid of an 
FEM connectivity scheme. Some results for a circular stacked-patch antenna 
will be presented to show the ability of the technique, and typical computed 
3-D near field distributions around the patches inside the two dielectric layers 
are displayed for visualization which reveals the insight of the EM coupling 
mechanism between the two patches. 


II. Hybrid FEM Formulation 

Consider the configuration of a stacked circular patch antenna shown in Fig- 
ure 1, where a cavity encompassing the patches is recessed in the ground plane 
to simulate the truncation of the computational region V. As described in 
(6), the finite element method (FEM) is employed to model the fields in the 
three dimensional volume, whereas a boundary integral equation on the aper- 
ture of the cavity (excluding the metallic patch area) is built up based upon 
the equivalence principle to associate the radiated fields in the region above 
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the ground plane to those in the cavity so that the continuity conditions of 
the tangential fields across the aperture is inherently imposed. (Interested 
Readers are referred to [6] for detailed discussions and mathematical deriva- 
tions.) It should be noted that the edge- based tetrahedral elements were 
used for the 3-D cavity fields and the triangular surface elements for the 
aperture fields simulation. The reasons are twofold (a) the edge elements 
have shown a precise representation of the vectorial electromagnetic fields 
without spurious modes; (b). arbitrarily shaped geometries and more so- 
phisticated feeding schemes (e.g. aperture coupled stripline, microstrip line, 
equi-potential model for coaxial cable (8], etc.) are easier to modei using 
edge elements. After the standard discretization and assembly procedures, 
we obtain a system of equations in a matrix form as 

N v N s 

E M {«?} +EI%) {9} -{**•>. t-1.2 N v (1) 

i j 

where i v and i‘ represent the unknown indices in the cavity and on the 
aperture, respectively, and the former contains the latter in global numbering 
since the aperture fields E* are part of the cavity fields E”. Nv and Ni are 
the number of cavity and aperture edges that do not lie on any PEC surfaces. 
It is apparent that the total number of the system unknowns is Nv and, in 
most cases, Ns <?C Nv- This property is important because [A,«j] is a sparse 
system although is a full matrix. Thus the storage requirement is 

insignificant for this exact boundary truncation, which governs the modeling 
accuracy for the radiation problems. 

In the case of stacked patch antennas, the lower patch is usually placed 
close to the cavity’s base, which is treated as a grounded PEC wall. He^ce, 
inside the FEM cavity, the lower patch should also be modeled as a PEC 
plate. Two separate dielectric layers above and below this patch must be 
considered to account for the substrate effects. Indeed, an attractive fea- 
ture of the FEM over integral equation method is its capability to model 
multilayer structures with the same case as it does with *.he homogeneous 
structures. 


III. Feed Simulations 

Probe Feed 

For the radiation analysis of a probe-fed patch antenna, we use the infinitesi- 
mally thin current source which is chosen to coincide with some mesh element 
edge. In this fashion, the excitation term {A',.] is usually modeled such that 
all the entries are set to zero except for those coincident with the probe(s) 
which link the patch and the ground plane. Specifically, if the probes are 
assumed to have a constant current /,,i = 1,2, M p (M p is the number of the 
probes), then the excitation entries associated with the corresponding edges 
may be written as 

K, = - jkZolili , ( 2 ) 

where k is the free space wave number, Z 0 is the free space intrinsic wave 
impedance and /, is the length of the i th edge. The phase difference between 
different probes may be readily represented by the phasor notation of each 
complex current /,-. Once the system of equations are solved, the input 
impedance is found to be 



where denotes the calculated field along the i th edge of length The 
radj'tion pattern may be obtained by integrating the computed aperture 
fields with the far— zone approximation of the dyadic Green’s function in the 
presence of an infinitely large PEC ground plane. 

Microstrip line Feed 

When the stacked patch antenna is led with a microstrip line network under- 
neath the ground plane (cavity’s base) via a coupling aperture (see Figure 2), 
special treatment for the feed structure must be considered in the FEM for- 
mulation. This is because the microstrip line is usually designed to have 
different size and shape as compared to the cavity’s geometric parameters. 
In addition, the microstrip line used as a transmission line, in most cases, 
carries traveling waves, thus, changing the field behavior in the cavity. Hence, 
the conventional FEM simulation treating the entire 3-D domain in a same 
fashion is not appropriate for this problem. 
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Consider the structure to be modeled as shown in Figure 3. In this case, 
it is appropriate to separate the computational domains because of the small 
element size required in modeling the feed structure. One difficulty encoun- 
tered when this decomposition is implemented is how one must model the 
coupling through the aperture. (We consider in this report a rectangular 
aperture which has been extensively employed in practice. Other aperture 
shapes for wideband consideration have also been seen in the literature and 
this formulation can be readily extended for those designs.) As described in 
[6], the cavity fields are discretized using tetrahedral elements, whereas in 
the microstrip line rectangular brick elements are considered the best candi- 
dates for field modeling. The selection of brick elements is particularly useful 
because of the geometry and the field property of the microstrip line, and it 
is shown that this choice retains better accuracy and efficiency in comparison 
to other element selections. As shown in Figure 4, although both types of 
elements employ edge-based field expansions, the meshes across the common 
area (coupling aperture) are different, and this causes the inconsistency for 
the field representation in the formulation. In the implementation, this im- 
plies that the usual boundary condition for the tangential field continuity 
across the aperture can not be enforced in a straightforward manner. 

We observe that the aperture is narrow in shape along the microstrip line 
and this fact allows us to assume a ‘static’ field feature for a given frequency. 
Therefore, the ‘potential’ concept may be applied to the problem in order to 
relate the fields below and above the aperture. Specifically, we enforce the 
‘equi-potential’ continuity condition, which states that: the potential across 
the narrow aperture due to the field distribution in the cavity region must be 

equal to that due to the microstrip line excitation. If the edges are classified 
as 


in mesh (b): 

• Ej 1 j — 1,2,3,... vertical edges 

• Ef j = 1,2,3,... diagonal edges 

in mesh (c): 

• Ej j — 1,2,3,... vertical edgesonly 

then the equi-potential’ continuity condition is expressed as 

E i' = ‘i E j (4) 
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(5) 


E f = ^(cyEJ + ty tl EJ +1 ), 

where 



t and d are the lengths of vertical and diagonal edges, respectively. (Thus, t 
happens to be the width of the narrow rectangular aperture.) 

The connectivity scheme for entirely different computational domains 
may be extended by generalizing this concept. It is apparent that this ap- 
proach makes the FEM implementation straightforward for different geom- 
etry/size domains that would be significantly inefficient if only one type of 
elements were used for modeling the structure. In addition, this technique 
ensures the system condition since the number of distorted elements in the 
mesh are minimized. 


Results 

As an example, a circular stacked patch antenna was used as shown in Fig- 
ure 1, where the geometry parameters were designed for the antenna to be 
operated in the frequency range from 7 to 9 GHz. When excited using a 
microstrip line, it is observed that the field across the coupling aperture has 
a Gaussian-like distribution, as expected. Figure 5 shows the gain com- 
puted as a function of frequency and as seen, the design is indeed wideband, 
as promised. Figure 6-9 show the radiation patterns for both polarizations 
in the E-plane and H-plane, respectively, and Figure 10 plots the input 
impedance, whose VSWR is 2:1, corresponding to a bandwidth of 15%. In 
Figure 11, we show the magnitude of the computed field distribution on the 
antenna’s radiating aperture when a probe feed was employed for excitation. 
It is noted that the energy is concentrated on the two opposite sides of the 
top patch and that the PEC cavity wall has little effect on the field near 
the antenna edges. Also, we observe from Figure 12 that the field couples to 
the lower patch from two locations even though only one probe is used for 
excitation. 
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Figure 1: Illustration of A Cavity-backed Stacked Circular Patch Antenna, 
which could be fed by a probe or by an aperture coupled microstrip line. 




FEM Mesh Terminated by 



6764 cm 



rj- 0.65 cm slotwidth«0.04cm 

r 2 = 0.75cm Ig- guided wavelength 

t- 0.08128 cm 

d-t- 0.18128 cm 

d<> 0.26256 cm 

s- 0.2 145cm 

D a 2.6 cm* * Calculations were done using this cavity size. 

Thus, performance will remain unchanged 
when army elements are separated by comparable 
distances. Tighter placement requires further 
investigation to maintain performance within specs. 

Figure 2: A stacked circular patch antenna with a microstrip line feed. The 
design of geometric and dielectric parameters will be used for modeling. 
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Figure 3: Cross-section of an aperture coupled patch antenna, showing the 

cavity region 1 and the microstrip line region II for two different FEM com- 
putation domains. 







Gain VS Frequency for a Stacked Circular Patch Anttena 




frequency (GHz) 


Figure 5: The gain pattern vs. freqency in GHz for circular stacked patch 
antenna. 
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lagnitude distribution on the antenna’s radiating aper- 
the energy is concentrated on two opposite sides of the 
ie PEC cavity wall has little effect on the field. 









